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As of early 2009, latest results on Galactic sources (mainly shell-type and plerionic supernova 
remnants), as observed in the very-high-energy 7-ray domain, are reviewed. A particular 
attention is given to those obtained with the H.E.S.S. experiment during its Galactic Plane 
Survey which now covers the inner part of the Milky Way. From the well identified 7-ray 
sources to those without any obvious counterpart and the putative Galactic diffuse emission, 
this observational window fully deserves to be celebrated during this International Year of 
Astronomy, as a new mean to image the Galaxy and reveal sites of particle acceleration, 
potentially at the origin of Galactic cosmic rays. 

1 Introduction 

Current generation of Imaging Atmospheric Cherenkov Telescopes (hereafter, lACTs) have re- 
cently revealed a new population of Galactic sources emitting in the very-high energy (VHE; E 
> 100 GeV) gamma-ray domain^^Sl^ j,^ particular, the H.E.S.S. experiment, through a Galactic 
Plane Survey performed over the last five years and covering the inner Galaxy {I E [-90°, 60°], 
|b| < 3°, see Figure [TJ, has accomplished a major breakthrough by revealing most of these new 
VHE 7-ray sources, as shown in Figure [2j A variety of source classes, identified (i.e. coincident) 
with sources known at traditional wavelengths, was found, among them several shell-type super- 
nova remnants (hereafter, SNRs), isolated or interacting with the surrounding medium, many 
young and middle-aged offset pulsar wind nebulae (hereafter, PWNe), some young massive star 
clusters, and a bunch of 7-ray binaries. In regards with SNRs and PWNe, one of the main 
pending question concerns the nature of the observed VHE 7-ray emission, which relates to the 
parent population of accelerated particles, or, in other words, the difficulty in disentangling the 
hadronic and leptonic contributions to the observed emission. This is in turn intimately linked 
to the more general question of the origin of Galactic cosmic rays (hereafter, CRs). As we shall 
see in the following, these questions can be efficiently addressed through a detailed investigation 
of the broadband spectrum of these sources, from radio to VHE 7-ray domains, coupled with 
the recent theoretical developments of acceleration mechanisms. 




Figure 1: H.E.S.S. significance image of the inner part of the Galaxy (£ £ [-60°,40°], |b| < 3°), as of 2008 (from 
Chaves et al., H.E.S.S. collaboration, 2008). The color transition from blue to red is set to 5 a. 



Besides these sources whose nature is firmly estabhshed thanks to the existing multi-wavelength 
observations, many others fall into the category of the so-called dark sources {i.e. with no clear 
counterpart at other wavelengths) This can be first explained by the fact that the majority 
of the VHE 7-ray sources are extended, on scales of the order of tens of arcminutes, with no 
clear sub-structure. Although current lACTs have reached unprecedented sensitivities and an- 
gular resolutions, the morphology of most of the faint sources can not be characterised precisely. 
Moreover, instruments in other domains (radio, infrared. X-rays) usually feature angular resolu- 
tions at the arcsecond / subarcminute scales, often coupled with relatively small field of views, 
which (1) does not permit one to perform deep surveys of the whole Galactic Plane, and (2) 
makes it difficult to reveal large-scale structures coincident with these newly discovered sources. 
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Figure 2: A history of VHE Galactic astronomy. The number of VHE Galactic sources is shown against the 
date, which corresponds to either the date of the publication or the date of the conference where the discovery 
has been officially announced. The number of Galactic sources has tremendously increased by steps over the 
last five years, particularly through the Galactic Plane Surveys (GPS, EGPS standing for the extension of the 
GPS) conducted with H.E.S.S.. HD-Gamma08 corresponds to the 4"" International Meeting on High Energy 
Gamma-Ray Astronomy which was held in Heidelberg, July 2008. The sources marked in red have been revealed 
in between these steps, mainly thanks to dedicated observations. The extrapolation in time depicted by the 
dashed line serves as an estimate of the number of new sources which might be revealed in the incoming era of 
H.E.S.S. n (Southern Hemisphere), MAGIC II and VERITAS (Northern Hemisphere). 



In this contribution, latest results on VHE Galactic sources are reviewed, with a particular 
attention to those obtained with H.E.S.S.. The well-identified cases, such as shell-type SNRs 
and PWNe, are discussed in sections [2] and [3] respectively, together with the implications and 
new questions related to the acceleration mechanisms and the nature of the VHE 7-ray emission. 
Some interesting cases of dark sources are exposed in section U and the origin of the putative 
VHE Galactic diffuse emission is discussed in section O The conclusion focuses on the interest 
of population studies, through the log N - log S distribution of all the VHE Galactic sources 
known so far, and the perspectives with the next generation of lACTs (cta1S31, AGislSSl). 



2 Shell-type supernova remnants 

Since the first speculation of Baade & Zwicky in 1934, the question of SNRs as the main sources 
of Galactic CRs up to the knee (~ 3 PeV) and beyond is not yet settled, in spite of several decades 



of important observational^^ and theoretical investigations^^^^^^. The broadband spectrum of 
these sources, from radio to VHE 7-ray domains, is the signature of particles accelerated at 
the shock fronts and radiating photons through several channels (synchrotron SC, non-thermal 
bremsstralhung, inverse-Compton IC and Tr'^-decay). Therefore, it represents by far our best 
access to the acceleration processes in 

SNR slSSl. Up to now, five pre yiously kno wn shell-t ype 
SNRs, namely Cas aSEI, RX J1713.7-3946[I2I29I35]^ J0852-4622[I3IEE3l36l^ rcW 86^, 
and more recently SN lOOG^SSl, have been discovered in VHE 7-rays, the last four exhibiting a 
shell-type morphology matching that observed in X-rays. As discussed before, the main question 
concerns the nature of the VHE 7-ray emission from these shell-type SNRs (leptonic through IC 
emissi on or hadronic thr ough vr'^-decay decay from p-p interactions) and gives rise to an intense 
debate E^^^'^ l ^^^^^ ^ ^^ l ^^i On one hand, the correlation between X-ray and VHE 7-ray emission 
would support the leptonic scenario but implies, in a simple one-zone model, a spatially-averaged 
magnetic field of the order of a few tens of fiG (except Cas A). This value seems uncomfortably 
low in comparison with the theoretical prediction of magnetic field amplication associated with 
the efficient production of CRs at forward shocks, by the (non-linear) diffusive shock acceleration 
(DSA) mechanism l^l^^l. Magnetic fields of > 100 have been derived from the measured 
thickness of the X-ray filaments in several youn and, more recently, from the fast 

jooi 

variability of small-scale X-ray dots and clumps'^, in case these localized structures effectively 
reflect the SC losses of high-energy electrons in strong (amplified) magnetic fields 0. On the other 
hand, for the four resolved shell-type SNRs, the lack of clear correlation between the tracers of 
the ISM and the VHE 7-ray emission, together with the tight constr aints on the local density 
derived from the absence (or the faint level) of thermal X-ray emission^EH], does not permit one 
to draw firm conclusions in favor of an hadronic origin in these shell-type SNRs§. 



Table 1: Observational constraints on VHE shell-type SNRs. Distances dsNR assumed here are 1, 1, 1, 2.2, 3.4, 4.8 
and 2.3 kpc, for RX J1713, Vela Jr, RCW 86, SN 1006, Cas A, Kepler and Tycho, respectively. The first column 
gives the magnetic field values assuming a one-zone leptonic model, with standard seed radiation fields (CMB and 
Galactic infra-red and star-light emissions) . The second column shows the magnetic field values derived from the 
thickness of the X-ray filaments. Third and fourth colums give the widths of VHE shells and of X-ray filaments, 
respectively, t^cr in the sixth column represents the fraction of the energy of the explosion, E51 (in units of lO^'^ 
erg), injected into CR protons, at the distance dsNR and for a gas density Ua = n/a cm"-^. Such density can then 
be compared to those given in the last column, mainly derived from the level of thermal X-ray emission. 
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Table [J summarizes the relevant parameters of the five shell-type SNRs detected in VHE 
"However, it was proposed that the thickness of the X-ray filaments would trace the magnetic fiel d d amping 
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downstream of the forward shock and, therefore, would not be a measure of the magnetic field strength' — As for 
dots and clumps, which would refiect the inherent turbulence of the magnetic field, even in the case of a steady 
particle distribution^^. 

''Note that, in the framework of the non-linear DSA, the post-shock gas temperature, which is expected to 
lie in the X-ray domain for young SNRs, and consequently the the rmal X-ray emission, can be reduced (shifted 
towards lower temperatures) due to strong particle acceleratio n^^ . However, there are still discrepancies on the 
level at which the thermal X-ray emission could be suppressed'^. 



7-rays and of the t wo historical SNRs, Kepler and Tycho, for which upper limits have been 
obtained so farEIMl. It appears clearly that the magnetic field values estimated from a one- zone 
leptonic model (fx/fvHE oc B^), are significantly below those derived from the thickness of the X- 
ray filaments. The situation looks even worse in the case of Kepler and Tycho SNRs, for which 
only rather strongly amplified magnetic fields seem to be compatible with the non-detection 
of VHE 7-rays, within these simple assumptions'^. On the other hand, the typical width of 
the X-ray filaments, over which amplified magnetic fields of about hundreds of //G may exist, 
are an order of magnitude below the widths of the resolved VHE 7-ray shells. Hence, in case 
the magnetic field has been damped quickly behind the forward shock, the observed emission 
could still be explained by IC emission, with a fairly weak spatially-averaged value. Moreover, a 
detailed modeling of the interstellar radiation field for the calculation of the IC spectrum may 
help to improve the fit to the VHE 7-ray spectra'^H]^ jj^ regards with the hadronic scenario, the 
energy injected into CR protons, required to explain the VHE 7-ray flux, is quite demanding, 
especially if the constraints on the gas density from the level of thermal X-ray emission are 
taken at face value. Apart from these energetical considerations, Tr'^-decay from p-p interactions 
seems to better explain the highest energy (> 10 TeV) data points measured in RX J1713.7- 
3946, where the Klein-Nishina limit of IC scattering takes place'^. Therefore, both scenarios, 
in their simpliest form, suffer from severe limitations. The question of these shell-type SNRs 
as efficient Galactic CR accelerators can only be efficiently addressed through spectro-imaging 
analysis in X-rays and VHE 7-rays, two domains whose current instrumental characteristics are 
quite different |, together with theoretical developments of the DSA mechanism. 



3 Plerionic supernova remnants 

Besides shell-type SNRs, a significa nt frac tion of VHE 7-ray sources is (or at least seems to 
be) associated with energetic pulsars ^-*- ' ^^ ! These sources can generate bubbles of relativistic 
particles and magnetic field when their ultra-relativistic wind interacts with the surrounding 
medium (SNR or interstellar medium) ^^Q] Their confinement leads to the formation of strong 
shocks, which can accelerate particles up to hundreds of TeV and beyond, thus generating 
luminous nebulae seen across the entire electromagnetic spectrum in the SC emission from radio 
to hard X-rays, and through IC process and potentially vr'^-decay from p-p interactions'^, in the 
VHE 7-ray domain. On one hand, recent advances in the study of PWNe have been made from 
mainly radio and X-r av o bservations of the complex morphology of the inner PWN structure 
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at the arcsecond scales'^. On the other hand, in the VHE domain, H.E.S.S. has proven to be 
capable to measure, in at least one case'SHj spatially resolved spectra at the tenths of degree 
scales. These complementary VHE observations then permit one to probe the electron spectra 
in these sources and to investigate the associated magnetic field distribution^^. 

Two classes of VHE 7-ray PWNe have recently emerged, based on observational grounds: 
young systems such as the Crab 

nebulalSni, GO.Q+O.lllS]^ MSH 15-52113 and the newly discovered 
VHE 7-ray sources associated with the Crab-like pulsars of 021.5-0.9^^ and Kes 75^^, and 
evolved (extended and resolved) systems (i.e. with characteristic ages Tc > 10^ yr), as exemplified 
by Vela X^, HESS J1825-137EI], HESS J1718-385 and HESS JlSOg-igS^I. in the former case, 
the VHE 7-ray emission, when resolved, matches quite well the morphology seen in X-rays, 
while in the latter case, these VHE 7-ray PWNe were found to be significantly offset from the 
pulsar position, with large size ratios between the X-ray and VHE 7-ray emission regions. The 
evolution of the SNR blastwave into an inhomogeneous and/or the high velocity of the 

pulsar's], together with a low magnetic field value (~ 5 fiG^^, may explain these large offset 



"^While soft X-ray (< 10 keV) telescopes feature angular resolutions at the arcsecond scales, both the imaging 
instruments above 10 keV and current lACTs reach angular resolutions of the order 5-10' at best. 



filled-center VHE 7-ray sources as being the relic nebulae from the past history of the pulsar wind 
inside its host SNR. Since VHE-emitting electrons are usually less energetic than X-ray-emitting 
ones, they do not suffer from severe radiative losses and the majority of them may survive from 
(and hence probe) early epochs of the PWN evolution. This interpretation has been further 
supported by the discovery of the spectral softening of the VHE nebula HESS J1825-137 as 
a function of the distance from the pulsarEB. Given the discovery of this large population of 
middle-aged many new sources regularly revealed by the on-going H.E.S.S. Galactic 

Plane Survey could fall into this category, some of them being classified as PWN candidates. 
For instance, HESS J1356-645^^1ies close to the young (tc = 7.3 kyr) and energetic (E = 3.1 x 
-|^q36 gj.g g-i) 166 ms pulsar PSR J1357-6429E!, for which only a marginal evidence of a 3" diffuse 
X-ray emission {i.e. of a putative PWN) was Interestingly, the extended VHE 7-ray 

source coincides with a diffuse radio emission, originally catalogued SNR candidate^. The 
on-going analysis of archival radio and X-ray data should thus help to constraint the nature 
of HESS J1356-645 (Aharonian et al., in prep.). Therefore, these sources can be confirmed as 
VHE 7-ray PWNe thanks to a detailed investigation of all the available multi- wavelength data, 
together with follow-up observations in radio {e.g. HESS J1857+026 - PSR J1856-h0245)H3Ml 
and with now Fermi/LAT {e.g. MGRO J1908-h06 - HESS J1908+063 - OFGL J1907.5-h0617) 
| 39 | 81 | 80 ]^ jj-^ order to reveal the associated (presumably energetic) pulsars. 

4 VHE 7-ray "nebulae" 

From the well-identified cases such as shell-type SNRs and PWNe, to the PWN candidates for 
which further data are required to firmly establish the putative association, it is presented here 
some cases of VHE 7-ray nebulae, also called dark sources^S. HESS J1731-347 represents the 
best example in this regard. Originally classified as a dark source, a faint and extended (R 
~ 0.25°) non-thermal radio and X-ray shell-type SNR, named G353.6-0.7, coincident with the 
extended VHE 7-ray emission, was discovered in the archival data'^. At a distance of ~ 3.2 kpc, 
estimated from HI absorption measurements toward an adjacent HII region, this SNR would 
have a physical diameter of ~ 28 pc, significantly larger than the known VHE shell-type SNRs 
described in section [2j This would then suggest that G353.6-0.7 is an old (~ 2.7 X 10^ yrs) 
and intrinsically very bright SNR. However, its distance and the nature of the X-ray emission 
are poorly constrained. Even though some theoretical studies have proposed that old SNRs (~ 
10^~^ yrs) could still emit in the VHE 7-ray domain^^, it is commonly thought that multi- 
TeV particles usually leave the acceleration site on timescales of a few thousands of years ^1^. 
Therefore, follow-up radio and X-ray observations are needed in order to shed further light on 
the nature of this newly discovered shell-type SNR. 

This example serves as a discussion about the VHE 7-ray emission from SNRs. During the 
Sedov phase of the SNR evolution, accelerated particles are gradually injected in the ISM, the 
most energetic ones being released first. In case the SNR lies close to a molecular cloud (MC, at 
< 100 pc), delayed VHE 7-ray (and neutrino) emission of the latter, through p-p interactions, 
may arise at detectable levels with the current 

lACTsSSl. 

The duration of VHE 7-ray emission 
from the cloud (> 10^ yrs) would then last much longer than that of the SNR itself, since 
it is determined by the time of propagation of CRs from the accelerator to the target. On a 
theoretical side, the detection of such emission would then indicate that the nearby SNR was 
acting in the past as an effective Galactic CR accelerator or PeVatron @. On an observational 
side, in this recently revisited scenariJ^, (some of) the unidentified VHE 7-ray sources could be 



"^Thus, the cutoff measured in RX J1713. 7-3946 at ~ 20 TeV'^, which translates into an Emax of particles at 
> 100 TeV (well below the knee in the CR spectrum observed at Earth), would imply that RX J1713. 7-3946 was 
a Pe Vatron in the past and that the highest energy CRs have already been released in the surrounding medium 



indirectly associated with old SNRs, the 7-ray emission being produced during the interaction 
of escaping CRs with nearby MCs. One would then expect a correlation between the VHE 
7-ray emission and the tracers of molecular material (^^CO, ^^CO and masers in case the SNR 
shock encounters the MC), as it might be the case for the VHE 7-ray emission detected by 
H.E.S.S. toward the old SN Rs W4 1 (HESS J1834-08712H]) and W28 (HESS J1800-240/J1801-233 
S51), the CTB 37 complex EOESl^ and HESS J1745-303H61. However, it is worth noting that 
some of these nebulae could be instead VHE 7-ray PWNe^^, as discussed in section O the 
large lifetime of VHE 7-ray emitting electrons in low magnetic field environments (~ 20 B^^q 

— 1/2 

E^ TcV ^y^) makes the ratio of the VHE 7-ray luminosity to the X-ray luminosity an increasing 

function of the source age and size^^. Therefore, one would expect VHE 7-ray PWNe to be 
hardly detectable with current X-ray instruments, and more generally at any other wavelength 
(leading to a VHE-on/?/, dark, source). Moreover, for most of these nebulae, the morphology 
is poorly characterized, where only the source barycenter and gaussian extension are usually 
provided. Many of them may well be multiple sources of different kinds, as for HESS J1800-240 
(A, B, & C)S5land HESS J1745-303SS]. Observations with the next generation of lACTs such 
as CTA and AGIS, with better sensitivities and angular resolutions, will undoubtedly help to 
search for counterparts with small field-of-view instruments and, thus, constrain the nature of 
the VHE 7-ray emission (s). 

5 VHE Galactic diffuse emission... or unresolved sources? 

This section is devoted to the recent detection by Milagro of a large-scale VHE 7-ray diffuse 
emission (30° < i < 110° and 136° < i < 216°, |b| < 10°)EI1, after removing the contribution 
of the sources detected by this experiment'^Sl. in the inner part of the Galaxy (between 30 and 
65° in longitude, i.e. excluding the peculiar Cygnus region), the traditional GALPROP model 
fails to fit the measured flux at ^ 15 TeV. An optimized version of GALPROP has been designed 
to reproduce the EGRET data by relaxing the restriction of the local CR measurements. In 
this model, the electron spectrum is contrained by the EGRET data themselves, such that any 
hard and faint (relatively to the standard vr'^-decay spectrum from CR protons) IC spectrum, 
as proposed by the authors, would not violate the GeV measurements, while explaining the 
measured flux at 15 TeV. Diffuse emission would thus be almost entirely explained by ~ 100 TeV 
electrons scattering off the CMB, with a flux, after propagation, of four times the one measured 
locally. Interestingly, this region has also been surveyed by H.E.S.S., featuring much better 
sensitivity and angular resolution than Milagro at the expense of a much smaller field of view 
(see Figure [3]), though with a non-uniform coverage'^. So far, five VHE 7-ray sources have 
been detected by H.E.S.S. (and unresolved by Milagro) in this region, and the resulting summed 
spectrum is shown in Figure [3l together with the Milagro measurement. Roughly 20 % of the 
diffuse emission is already explained by these H.E.S.S. sources, and a larger fraction should be 
reached in a next future once the existing H.E.S.S. data will be carefully analyzed and the survey 
will become uniform. First studies of the VHE 7-ray source population, based on the second 
H.E.S.S. survey catalogue^, had already suggested that at least 10 % of the VHE Galactic 
diffuse emission should be attributed to unresolved sources^^. 



6 Conclusion 

To conclude this review, the log N(>S) - log S distribution of all the VHE Galactic sources known 
so far is shown in FigureHl A slope of -1 corresponds to a uniform infinite plane distribution while 
a slope of ~ -0.6 indicates that the population follows the spiral arm distribution on Galactic 
scales {i.e. at distances of > 8kpc), as expected from young stellar population generating SNRs 
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Figure 3: Flux - Size plot of the VHE Galactic sources detected so far. The integrated flux is given between 1 and 
10 TeV, in units of the Crab nebula in the same energy band. The size correponds to either the intrinsic mean 
source width or, in the case of shell-type SNRs, the outer radius, in units of degrees. The red points represent 
the sources detected by Milagro during its Galactic Plane survey {£ £ [30,65]°, |b| < 2°), while the blue points 
correspond to those revealed by H.E.S.S. in the same region of the sky. The solid lines show the sensitivities 
of the two instruments, according to their respective characteristics given in brackets, which degrades with the 
source extent as: S = So x (R^ + a^^^)^^^ /ap^f, where So is the nominal sensitivity and Rs is the effective source 
size. Note that this law does not hold anymore when the source size becomes comparable to the instrument field 
of view, which explains why the H.E.S.S. sensitivity curve is valid only for Rs < 1°. The inset plot shows the 
Milagro diffuse emission measured at 15 TeV, the expected diffuse flux from the optimized GALPROP model and 
the summed spectrum of all the H.E.S.S. sources falling into the region probed by Milagro. 



and PWNe. Even though there are inherent comphcations in the interpretation of such plot, the 
transition between the two regimes (at ~ 0.15 Crab) seems to take place above the completeness 
limit of the H.E.S.S. survey. Extrapolating the curve down to 1 mCrab in sensitivity then 
implies that the next generation of lACTs like CTA^SH ^nd AGIS^^ should detect at least 
500 sources. As discussed in section HI many faint sources could actually be several sources, 
not clearly resolved yet, leading the log N(>S) - log S to soften toward low fluxes. The paper 
focused mainly on the latest results and open questions related to shell-type SNRs and PWNe. 
Besides them, new classes of VHE 7-ray emitters in the Milky- Way are expected to emerge, 
either from multi-wavelength follow-up observations of the so-called dark sources, as it might be 
the case of HESS J1503-582ES]^ or through dedicated observations of sky regions of interest, as 
exemplified by the recent detection of VHE 7-rays towards the young stellar cluster Westerlund 2 
^33. Such detection has triggered many exciting questions on the nature of the VHE 7-ray 
emission and more generally on the contribution of such wind-blown bubbles to the Galactic 
CR flux. HESS J1848-018 could be the second case of this kind, where the VHE 7-ray emission 
is found to be slightly offset from the massive star- forming (Galactic "mini-starburst" ) region 
W 43^^. Moreover, 7-ray binaries have now joined the club of VHE 7-ray emitters and many 
pending issues still need to be investigated. Next generation of VHE observatories, in tight 
link with incoming multi- wavelength instruments, will undoubtedly play a crucial role in the 
understanding of all these acceleration sites in the Galaxy and, very likely, in the discovery of 
an even larger diversity of VHE 7-ray sources than expected. 
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Figure 4: Log N(>S) - log S diagram of the VHE Galactic sources. The dashed lines represent the dispersion of 
the distribution after taking into account the statistical and systematical errors on the source spectra through a 
Monte-Carlo simulation. The red and orange solid lines correspond to the expectations from a uniform distribution 
in a thin disk (slope of -1) and a source population distributed along the spiral arms (slope of -0.6), respectively. 
The transition between these two regimes was set arbitrarly. The dot-dashed line shows the completeness limit 
of the H.E.S.S. survey (for a source extent of 0.2° in radius, and five hours of observations everywhere in the 
inner part of the Galaxy). The extrapolation down to 1 mCrab, a sensitivity that should be reached by the next 
generation of lACTs like CTA and AGIS, would then lead to the detection of about 500 VHE Galactic sources. 
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